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Abstract. In this article we study the asymptotic behavior of incompressible, ideal, time- 
dependent two dimensional flow in the exterior of a single smooth obstacle when the size of the 
obstacle becomes very small. Our main purpose is to identify the equation satisfied by the limit 
flow. We will see that the asymptotic behavior depends on 7, the circulation around the obstacle. 
For smooth flow around a single obstacle, 7 is a conserved quantity which is determined by the 
initial data. We will show that if 7 = 0, the limit flow satisfies the standard incompressible Euler 
equations in the full plane but, if 7 7^ 0, the limit equation acquires an additional forcing term. We 
treat this problem by first constructing a sequence of approximate solutions to the incompressible 
2D Euler equation in the full plane from the exact solutions obtained when solving the equation on 
the exterior of each obstacle and then passing to the limit on the weak formulation of the equation. 
We use an explicit treatment of the Green's function of the exterior domain based on conformal 
maps, a priori estimates obtained by carefully examining the limiting process and the Div-Curl 
Lemma, together with a standard weak convergence treatment of the nonlinearity for the passage 
to the limit. 
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1. Introduction 



Let f2 be a bounded, connected and simply connected domain of the plane with smooth 
boundary F and, for each e > 0, let fi^ = eQ with boundary F^. We consider a family 

= u^{x,t) of incompressible flows which satisfy the two-dimensional Euler equations on 
the exterior of fig, with velocity tangent to F^ and satisfying the initial conditions (i) the 
initial vorticity u'^^x) = curl u^(x,0) is independent of e and the support of does not 
intersect the origin and (ii) the circulation around F^ of the velocity is a real constant 7 
which does not depend on e either. Our purpose in this article is to identify the asymptotic 
behavior when e — > of the sequence {m^}. We will prove that for 7 = 0, converges 
to a flow u satisfying the incompressible two dimensional Euler equations in the full plane 
with initial vorticity and that, for 7 7^ 0, any weak limit of the sequence {u"^} satisfies 
a modification of the Euler equations which, in the vorticity formulation, takes the aspect 
of an additional convection term. 

In order to clarify the issues involved in this problem, let us consider that, for each 
fixed time, the velocity u of the flow around a small obstacle decomposed as u = Uo + Ub, 
where Ub is the background velocity, which is slowly varying with respect to the scale of 
the obstacle, incorporating what we imagine the flow velocity would be if the obstacle was 
not present, and Uo is the correction to Ub due to the presence of the obstacle. As the 
obstacle disappears towards a point p, the background flow Ub appears as a constant flow 
background Ub{p) relative to an observer on the obstacle. One can see that, when Ub{p) 7^ 0, 
the fact that the full velocity u is tangent to the obstacle means that Uo has to be a large 
perturbation in absolute value. On the other hand, one expects the perturbation produced 
by the small obstacle to be sharply localized, so that \uo\ should converge pointwise to zero 
away from the obstacle. The small obstacle generates large velocity gradients in a nearby 
region. How this effect influences the limiting process is the main point of the present 
work. 

It is a well known fact that the ideal flow assumption is physically inappropriate to 
model the behavior of the flow near an obstacle, due to boundary layer effects. Hence, the 
whole issue of small obstacle asymptotics would be more physically meaningful if posed 
for the Navier-Stokes equations. The asymptotics of the ideal flow case which we present 
here should be regarded as a first step towards the rigorous analysis of the small obstacle 
problem. 

Most of the work on time-dependent, incompressible exterior flow has been in the nature 
of well-posedness through energy methods, see P, |^ . Energy estimates are global, so that 
it would be difficult to treat our sharply localized problem through such means. The 
alternative is to adopt the vortex dynamics point of view. This means understanding 
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2D flow in terms of tlie description of the dynamics of vorticity, an approacli wliicli lias 
been very fruitful recently. The most important recent results on 2D Euler are Chemin's 
Theorem, on the regularity of vortex patches, see [0 , which is by nature a vortex dynamics 
result, and Delort's Existence Theorem, see 0], which relies heavily on the vortex dynamics 
point of view. 

Our interest in the small obstacle problem was motivated by the problem of confinement 
of vorticity. Let us examine briefly the nature of this problem and survey some of the 
results obtained thus far. Consider uj = uj{x,t) a (classical) solution of the full plane 2D 
Euler equations such that uj{x, 0) is compactly supported. The problem of confinement 
of vorticity is to obtain control over the growth of the diameter of the support of uj{-,t). 
Current research on this subject originates with a result obtained by C. Marchioro in ||TU|. 
He proved that the solution of the incompressible 2D Euler equations in the full plane, 
with bounded, nonnegative initial vorticity with support contained in the ball B{0;Rq), 
will have, at time t, its support contained in a ball of radius R{t) = {ht + for some 

constant 6 > 0. The state-of-the-art confinement result in the full plane has almost fourth 
root exponent, see 0. In [^, Marchioro addressed the problem of confinement for 
exterior flow. Using the techniques developed for full plane flow, Marchioro proved cubic- 
root confinement for the exterior of a disk and almost square root confinement for a general 
exterior domain. We will not get into the issue of why the confinement result is sensitive 
to the presence and to the geometry of the domain, but we observe that confinement is 
connected to the way an obstacle influences very distant particles, called far-field effects. 
The scaling behavior of the incompressible Euler equations makes the problem of describing 
the influence of a very distant obstacle for a very long time naturally associated to observing 
the effect of a vanishingly small nearby obstacle for some time, which is the object of this 
article. In this context, the influence of the precise shape of the small domain in the vortex 
motion is of particular interest. 

From the technical standpoint, this article makes use of the techniques of weak con- 
vergence methods for the asymptotic analysis. Such methods have often been developed 
for proving existence of weak solutions, see for example 0, H], but they are well suited 
for studying singular limits in general. The basic ingredients of the present proof are a 
collection of a priori estimates obtained mainly through exhaustive use of explicit formulas 
for the Green's function of the exterior domain and strong compactness of approximate 
velocities obtained by using the parametrized div-curl Theorem introduced in [Q. 

The remainder of this paper is divided into five sections. In the second section, included 
mostly for completeness' sake and for fixing notation, we collect information on classical 
potential theory for the Laplacian in an exterior domain and on holomorphic maps as 
required to write explicit formulas for the Green's function in terms of the Riemann map 
associated to the obstacle. In the third section we formulate precisely the exterior flow 
problem and the problem of small obstacle asymptotics. In the fourth section we derive the 
collection of a priori estimates required for the passage to the limit, identifying a collection 
of quantities that remain under control in the small obstacle asymptotics. We prove our 
main result in the fifth section, identifying the PDE in the full plane satisfied by the limit 
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flow. In the sixth and last section we collect our conclusions and point the way for future 
investigation on this problem. 

2. The Laplacian in an exterior domain 

The purpose of this section is to collect a number of facts associated to classical potential 
theory for the Laplacian on an exterior domain in the plane. We claim no originality 
on these results, although the explicit form in which they are presented is not the one 
commonly found in the literature. We include a thorough discussion here for the sake of 
completeness. 

2.1. Conformal maps. Let Qhe a. bounded, open, simply connected subset of the plane, 
whose boundary, denoted by F, is a C°° Jordan curve. We will denote by 11 the unbounded 
connected component of - T, so that Q"" = TT. Let D = B{0; 1), S = dD. 

In what follows we identify with the complex plane C. We begin by constructing a 
smooth biholomorphism between 11 and {a; G | \x\ > 1} = int D'^. 

Lemma 2.1. Let 11 = int VL'^ he as above. There exists a smooth biholomorphism T : 
n int D'^, extending smoothly up to the boundary, mapping F to S . Furthermore, there 
exists a nonzero real number (3 and a bounded holomorphic function h : U ^ C such that: 

(2.1) T{z) = pz + h{z). 
Additionally, 

(2.2) h'{z) = O 

Proof. We may assume without loss of generality that -8(0; 1) C Q. If this were not the 
case then suitable dilation and translation would make it so. ^ 

For z e n, consider the holomorphic map given hj z ^ I{z) = 1/z and denote by Q, 
F, the image of 11, F respectively, under this map. There exists a biholomorphism TZ from 
i7 U {0} to 5(0; 1) given by the Riemann Mapping Theorem, mapping F to S. Since F was 
assumed to be a C°° Jordan curve so is F, hence TZ and its inverse 7l~^ have continuous 
extensions up to the boundary, along with all their derivatives (see ||l|). We choose the 
Riemann mapping TZ such that 7^(0) = and we define T to be given by the compositions: 
T = I olZo I . It is easy to see that T is a biholomorphism, continuous up to F, along with 
all of its derivatives. 

Consider the function / defined by 

TZ{z) = zf{z), zeh. 

This is a holomorphic function on Q which cannot vanish at since 7^ is a bijection. By 
the same token, / cannot vanish anywhere else in the closure of fl. It follows that 1/ f{z) is 
also a holomorphic function. Hence we can write 1/ f{z) = (3 + zg{z), with g a holomorphic 
function and p = l//(0). We find g{z) = -{f{z) - /(0))//(0)7^(4, which can be seen to 
be bounded since 2; = is a zero of order 1 for TZ and / is bounded. 




4 



Thus we obtain 



T{z) = l/nil/z) = z{(3 + {l/z)g{l/z)) = (3z + g{l/z). 

The desired function h is given by h{z) = g{l/z). We may assume without loss of generahty 
that P is real by multiplying the function T constructed above by if necessary, which 
does not change the desired properties. Finally, we observe that 




and, by construction, g' is a holomorphic function in fl, bounded in the closure. □ 



Remark 2.1. It follows by construction of T that 

T'{z) = n'{i/z)/{f{i/z)r, 

which can be easily seen to be bounded from above and below. Therefore, if we make the 
canonical identification of with C, writing z = xi + 1x2 for the point x = {xi,X2), then, 
if DT stands for the Jacobian matrix associated to T', there exist a positive constant C 
such that 

(2.3) II^^IU- < C and HDT^^loo < C. 

2.2. Explicit formulas for the Green's function. Here we will obtain an explicit 
formula for the Green's function of the Laplacian in 11 in terms of the conformal mapping 
above. We denote this Green's function by Gn = Gu{x,y)] we must have AyGu_{x,y) = 
S{y — x), GYi{x,y) = for |/ G F and Gn{x,y) = Gu{y,x). If Q = B{0; 1) and if we write 
X* — x/\x\'^ then there is a unique Green's function which can be written explicitly as 

GDc{x,y) = 7^ log 

It is easy to verify that if Xq G (-8(0; l))"^ and h satisfies A/i = 5{x—Xq) in a neighborhood 
of Xq then h = hoT satisfies Ah = 5{y — T~^{xo)) in a neighborhood of T~^{xo). We use 
this fact on G to write: 

We will concern ourselves mostly with first order derivatives of the Green's function, 
which we will introduce through the notation Ku — Kii{x, y) = V^Gu{x, y). The explicit 
formula for is given by 

.2 .^ K ir v) = ((^X^O - Tiy))DT(,r))^ _ ((T(.r) - (l\y)r)DT{.r))^ 

'^^ 27r|r(x)-r(y)p 27r|r(x) - (r(y))*|2 • 

Note that Ku{x,y) = DT\x)KDc{T{x),T{y)). 





\x-y 




\x 


- y*\ 


\y\ 
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We require information on far-field behavior of Ku- We will use several times the fol- 
lowing general relation: 



(2.6) 



a b 



\a — b\ 



which can be readily checked by squaring both sides. 

We now find from the estimates (p.3|), (^TSl) and (p.6|), that 



( M<r \ny)-{ny)y 
\Kn{x,y)\ < C- 



\T{x)-ny)\\T{x)-iT{y)y\- 
For / G C^(n), we introduce the notation 

(2.7) Kn[f] = Kn[f]{x) ^ [ Kn{x,y)f{y)dy. 

Ju 

It is easy to see that the pointwise estimate for Ku above, yields, for large the 
estimate: 

(2.8) l^n[/]|(x)<^, 



where the constant C > depends on the size of the support of /. In the last inequality 
we have used the explicit formula for the biholomorphism T (|2.1| ). 



Lemma 2.2. The vector field u = K]j[f] is a solution to the elliptic system: 

div u = in U 
curl u = f inU 
u ■ n = on r 
limi^l^oo l^^l = 0. 

The proof of this Lemma is straightforward. 

2.3. Harmonic Vector Fields. The standard version of Hodge's Theorem is proved for 
compact manifolds, see [|16[ , with a natural extension to compact manifolds with boundary. 
Extending Hodge's Theory to noncompact manifolds is a difficult problem, see @] for a 
broad overview. There is a lot of special structure for the exterior domain, specifically the 
complex structure, which allows us to prove an elementary extension of Hodge's Theorem 
to this particular case. It is no surprise that this extension is valid, moreover, such a fact 
has been extensively used in the literature. As before, we include a complete discussion 
here because we will require the very explicit treatment involved for the remainder of this 
article. 

Let n be a bounded, open subset of whose boundary F is a smooth Jordan curve and 
let H = int fl'^. We will denote by n the unit normal exterior to Q'^ at F. In what follows all 
contour integrals are taken in the counter-clockwise sense, so that J^F -ds = — J^F ■ h^ds. 
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Proposition 2.1. There exists a unique classical solution H — Hu of the problem 

div H = 0, in 11, 
curl H — 0, in H, 

(2.9) I H ■n = Q onV, 

\H\ ^0 as X ^ oo, 

Moreover, Hu — 0{l/\x\) when \x\ — > oo. 

Proof. We will start by proving uniqueness. Let us assume there are two solutions Hi, H2 
and consider their difference H = Hi — H2- Since the problem is linear iJ is a harmonic 
vector field with zero circulation. We must show that H vanishes identically. Let us begin 
by estimating H at infinity. We will use the canonical identification of with C, writing 
z — xi-\- 1x2 instead ol x — {xi,X2). First note that H — 0{\x\~^) as x — > 00. To see 
this consider the holomorphic function / = f{z) — H{l/z). The singularity at 2; = is 
clearly removable and hence / is holomorphic in a neighborhood of the origin. Furthermore 
/(O) = 0, which implies that f{z)/z is also holomorphic in a neighborhood of 0. Thus we 
have shown the desired behavior at infinity. Next, note that, after an easy calculation we 
find, for any real vector field L, 

(2.10) jLdz = J L-ds-i J L-nds. 

Using this identity for H we conclude that the contour integral of ^ on F vanishes. Since 
H is holomorphic in Q'^, its contour integral on any closed curve enclosing Q vanishes as 
well. Let us consider the Taylor expansion of / = f{z) at ^ = 0, f{z) — cz + 0(^^). Fix 
s > suitably small. We have that: 

27ric = / ^^dz = / H{w)dw = 0. 

J\z\=e ^ J\w\=e^^ 

From this we conclude that actually H = 0{\x\^'^) as a; — 00. 

Next we observe that, since H is curl-free and has zero circulation, there exists a function 
(p such that H — V^p. This function is constructed by integrating along arbitrary paths. 
From the behavior of H at infinity and the construction of </? we have that \^{x) \ — 0{\x\''^) 
as a; ^ 00. 

Let us now use the estimates gathered above to show that H = 0. We integrate: 

/ \H\'^dx — lim / \H\^dx — lim / H ■ Vipdx — lim / (pH ■ —ds, 

Jn<= J B{0;R)\n J B{0;R)\n ^^°°J\x\=R R 

by integration by parts and using the fact that H is tangent to F. This limit vanishes 
because of the behavior of if and H at infinity. Thus we conclude that H is identically 
zero in Qf^. 



7 



We now obtain the existence of a classical solution. Let T be the biholomorphic mapping 
from Lemma |2.1| . Let ip = ip{z) = log |T(z)| for z G Q'^. Since the logarithm is a harmonic 
function in 5(0; l)'^ and T is analytic in LI it follows that ip is a. harmonic function of x. 
Furthermore ip = onT. Define U = V-^ip = {—dx^'^P-, dx^ip)- It can be easily checked that 
U satisfies all but one of the conditions in system ( |2.9| ), namely the condition on circulation 
on L. Since U is not identically zero it follows from the uniqueness part of the proof that 

[/ . ds = c ^ 0, 



r 

and hence we may take Hu = U/c. 

Finally we address the asymptotic behavior of Hu at infinity. Recall that, in the be- 
ginning of the proof of uniqueness above, we showed that the difference of two solutions 
behaved like 0(|x|~^) as x — > oo. The argument we gave did not depend on the circulation 
of the difference H. Therefore the same argument can be used with Hi = Hjj and H2 = 
to show the desired behavior for Hu- □ 

By construction the harmonic vector field Hn above is given by 

Hn{x) = CV^\og\T{x)l 

for some constant C. We argue that C = l/2n. Indeed, it is an easy calculation to prove 
that 

Huiz) = -Cti\ogTyiz) = -Ct^. 
Using (|2.1CI| ) together with the fact that Hu is tangent to L we find: 



Hu-ds= I Hudz = C I ^^dz = C I — = 2-kC, 
IT Jr Jr ^T[z) Js tw 

so that we must have C = 1/2tt in order to satisfy the last condition of system (|2.9| ). 
Throughout the remainder of this paper we set 

(2.11) Hu(.) = ^vHom.y i^^^w" i^r.(.)(rw)- 



2tt V 2n \T{x)\^ 2ti |T(x)|2 

3. Flow in an exterior domain 
The purpose of this section is to formulate precisely the small obstacle limit. 

3.L The initial-boundary value problem. We begin by formulating precisely the 
initial-boundary value problem for incompressible ideal fluid flow in an exterior domain. 
Let r be, as before, a smooth Jordan curve in the plane, dividing it into a bounded con- 
nected component, which we call VL and an unbounded connected component denoted 11. 
For x G r, denote by h{x) the exterior normal to VL'^ at x, as before. 

Let u = u{x,t) = {ui{xi,X2,t),U2{xi,X2,t)) be the velocity of an incompressible, ideal 
fluid in Q'^. We assume that u is tangent to F and m — when |x| 00. The evolution 



8 



of such a flow is governed by the Euler equations. We write the initial-boundary value 
problem as: 



(3.1) 



Ut + u- Vu = —Vp in n X (0, oo) 
div u = in n X [0, oo) 

u ■ h = on r X [0, oo) 

lim|a;|^oo \u\ = for t G [0, oo) 

in Q", 



u{x, 0) = Uo{x) 



where p = p{x, t) is the scalar pressure. If Uq is sufficiently smooth, global well-posedness 
of this problem was proved by K. Kikuchi in 0. 

Let LJ = curl u be the vorticity associated to this flow. In order to write a vorticity 
formulation of problem ( |3.1D we must be able to recover velocity from vorticity. The 
coupling of velocity and vorticity is given by the elliptic system 

div M = in n X [0, oo) 
curl u = u inU X [0, oo) 
u ■ fi = on r X [0, oo) 
lim|a;|^oo \u\ = for t G [0, oo). 

In view of Lemma and Proposition ^]T| this system has a unique solution up to a 
harmonic vector fleld, given by -u = u{x,t) = Ku[uj{-,t)]{x) + aHu{x), for some time- 
dependent function a = a(t). 

Lemma 3.1. If u is a smooth solution of ( \3. ]\ ) then a is constant in time. 
Proof. We introduce the stream function ip = il){x,t) given by: 



il){x,t) = Gu[i^]{x,t) = / Gu{x,y)uj{y,t)dy 



(3.2) 



Next, we observe that, by Kelvin's Circulation Theorem, the circulation of u around F 
is constant in time. Hence we have: 
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j u ds = j V-^ip ■ ds + a{t) j Hn-ds = a{t) - J Vijj ■ fids. 



Integrating by parts we flnd 
7 = a{t) 



lim 

R-*oo 



uodx 



B{0;R)\n 



Vip ■ ^ds 



dB{0;R) 



R 



a{t) — / iudx 



usmg 



Hence, since mass of vorticity is conserved, 
a(t) = 7 + / curl uodx 



□ 
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Finally, we can now write the vorticity formulation of this problem as: 

uJt + u ■ Vuj = 0, in n X (0, oo) 



(3.3) 



u = Ku[uj] + aHu, 
u{x, 0) = curl uo{x) 



in n X [0, oo) 
in n. 



3.2. The evanescent obstacle. In this subsection we will formulate a family of problems, 
parametrized by the size of the obstacle, in order to identify the asymptotic limit under 
consideration. Fix ujq = u)o{x) G C^(]R^) and assume that the origin does not belong to 
the support of ljq. Let be a domain in the plane satisfying the hypothesis of Lemma 
We will consider the family of rescaled domains fi^ = efl and we note that there exists 
such that, for all < e < Eq, the support of t^o does not intersect Q^- 

Let = dQ^ and = int Q'^. We denote the harmonic vector field given by Propo- 
sition ^Tl] by H'^ = Hu^. We also denote the Green's function for by and the 
corresponding kernel (and integral operator) K"^ = V"'"G^. 

Consider the system (E^) given below: 

+ ■ Vuj' = 0, in x (0, oo) 




^{x, 0) 

It follows from the work of Kikuchi M 



UJo[X) 



in X [0, oo) 
in n,. 



that, for any e > 0, if cuq is sufficiently smooth then 
this system has a unique smooth solution. 

In fact, Kikuchi's result asserts existence for the velocity formulation of (-Ee), which 
means that there exists a pressure such that u'^ and are a solution of problem ( |3.1| ) 
in the domain 11;. with initial velocity 

(3.4) ul = K'[uJo]+aH'. 



Consider T the biholomorphism from 11 to int constructed in Lemma |2.1| . Observe 
that T'{z) = T{z/e) defines a biholomorphism between 11^ and int V^, which extends 
smoothly up to the boundary and which maps F^ to S. We can use to write explicit 
formulas for K' and H', 



K' 



(3.5) 



and 



(3.6) 



X 



by recalling 

r 



and (pll) . We have 



-DT 



( mx/s)-T{y/s)))^ _ mx/e) - {T{y/e)y))- 
^ ' \ \T{x/e) -T{y/e)\^ l^r^ /.^ _ r^r,, /.^^*|2 



X 



2718 



DT\x/e) 



T{x/e) 

inx/e))- 
\T{x/eW 



iny/e)) 



We will require information on the behavior of H'^ as e — > 0. One easy observation on 
that regard is the fact that for any R> 0, there exists C = C{R) > such that 



(3.7) 



\L^{UenB{0;R)) 



<c, 
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uniformly in e. Indeed, 



\T{x/eW 

< [ -^^e'\det{DT-')iy)\dy 

Jl<\v\<C(R)/e ^ \y\ 



dx 



'l<\y\<CiR)/e ^ \y\ 

where we used the change of variables y = T{x/e), we used ( p.3|) and C{R) is a suitable 
constant computed using the expression for T from Lemma From this last expression, 
estimate (|3.7| ) follows easily. 

4. A PRIORI ESTIMATES 

In this section we will prove the a priori estimates on and tu'^ which are required to 
identify their asymptotic behavior. 

4.1. Velocity estimate. The key ingredients on the rigorous treatment of weak solutions 
for the incompressible Euler equations are usually the energy estimate on velocity and 
estimates on vorticity based on the fact that the vorticity is rearranged by incompressible 
flow, i.e. the distribution function of vorticity is a constant of motion. There are numerous 
instances of this observation, and we refer the reader to for a systematic discussion. For 
the problem in hand, the usual rearrangement estimates on vorticity hold for the sequence 
u^, due to the transport nature of equation {E^). For e small enough, so that the support 
of Uq is compactly contained on we have the a priori bounds 

(4.1) ||u;'^(-,t)||ioo(n^) = ||a;o||L°°(R2) and ||to'^(-, t) ||ii(n^) = ||ci;o||z,i(iR2) 

On the other hand, we are not able to use the energy estimate in the usual way. In order 
to explain what is at play here, let us illustrate the behavior of the velocity as £ — > with 
an explicit example. 

Example: Let = B{0;e), so that 11^ = {|a;| > e}. We consider the vorticity 

, , f 1 if 1 < |x| < 2 
^ ~ [ otherwise, 

which is a stationary solution of (E^) for any e < 1. 

The finite energy part of the associated velocity is = -ft'^[c<^], which, after straightfor- 
ward calculations gives 

Ifw if e < Ixl < 1 

2 \x\^ — I I 



v'^(x) 



if \x\ > 2 
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Clearly, there is (logarithmic) local blow-up of the kinetic energy as e ^ 0. The harmonic 
vector field associated to Ilg is the restriction of H{x) = x-^/27c\x\'^ to He, so that the full 
velocity is = + aH. If we want to single out a locally uniformly square integrable 
portion of for estimating we must add a well chosen part of aH to v^. We consider 

= (f^ + StcH) + {a — ?)Tt)H. Neither part of this decomposition is square-integrable at 
infinity, but v'^ + 37rif is not only uniformly bounded in Lf^^ but it is actually uniformly 
bounded in We observe that Svr = / uj{x)dx. The purpose of this subsection is to 
prove that something similar holds in general. 

Returning to the general sequence a;^, u^, assume again that e is small enough so that 
the support of c^o is compactly contained in lie- We introduce 



171= uo^dx = / Uodx. 

Introduce also 



re 



so that = f ^ + (a — m)H'^. 

Theorem 4.1. There exists a constant C > that depends only on the shape of Q such 
that 

II en ^ r^U £l|l/2 II £||l/2 

\\v \\l°°{u,) S <-^||t^ IIli • 

Proof. We write v'^{x,t) = Xi +X2 explicitly using (|3.5|) and ( p.6|) with 

27re V \T{x/e) -T[y/e)\^ 



and 



We begin by estimating Xi . Using ( |2.3| ) we get 



^ in, \T{x/e) -T{y/e)\ 

Let J = J(^) = I det(DT^^)(^)| and z = eT{x/e). Then, making the change of variables 
rj = eT{y/e), we find 

(4.2) |Xi| <C [ ^^\u;^{eT~\r]/6),t)\J{ri/e)dr]. 

J\v\>e \Z-V\ 



\>e \Z-V\ 

Next, we introduce 
(4.3) f%v,t) = \u;'{eT-\rj/e),t)\J{rj/e)x\rj\>e, 
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where xe is the characteristic function of the set E. We change variables back and we get: 

(4.4) ||r(-,t)|Ui(M^) = ||a;^(-,t)|Ui(M^), 
and, since by (|2.3| ) J is bounded, we find that 

(4.5) \\r{-ML^m<c\\u'\\L^. 

We deduce from (^) and (^) that 

1 



(4.6) 



< c 



This last inequality is Lemma 2.1 in According to ( |4.4|) and (|4.5|) this concludes the 
estimate for Xi. 

Let us now estimate X2. Using again ( ^731) we find 



IX, I < 



T{x/e) - {T{y/e)r 



T{x/e) 



\uj'{y,t)\dy. 



\T{x/e)-{T{y/e)y\^ \T{x/eW 

Again we consider J = J(^) = | det(DT^-'^)(^)| and z = eT{x/e) and we make the same 
change of variables r] = eT{y/e). Using ( pl6|) this yields 



IX, I < 



c 



< c 



z/e - er]* 



z/e 



\z/e — er7*P \z/e\'^ 
e^\r]*\ 



oj'{eT-\'n/e),t)\J{'n/e)d'n 



r,\>e \z\\z-e'^r]*\ 



iu'ieT-\7]/e),t)\Jir]/e)dr]. 



As 2; = eT{x/e) and the image of T is the exterior of the unit disk, it follows that |2;| > e. 
Hence, 

e\7]* 



IX2I < C 



r,\>e \z-e^7]*\ 



\uj%eT-\r]/e),t)\J{r]/e)dij. 



We again change variables in the integral above, writing 6 = erj*. Then we have: 

= c(f +[ )=C(In+I^i)- 

\J\e\<i/2 Ji/2<\e\<iJ 

First we estimate X21. For |^^| < 1/2 we have that I2; — e^| > so that 



|X2i|< / 2e\e\\u^{eT-\e*),t)\J{e*)^^ 
\e\<i/2 Fl 



r?|>2e 



\uj^{eT~\7^/e),t)\J{7^/e) 



dr]<2 



di], 
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where was introduced in ([4.3|) . Hence, it follows from ( |4.4| ), ( [4 .51) and ( [4.6|) with z = 

that IX21I < C||c<j^||^(^||co''^||]/c^. 
Finally, we estimate X22. Let 

g%9,t) = \co%eT-\9*),t)\J{9n^^ 

so that 

(4.7) J22= / -^^f{e,t)de. 

Jl/2<\9\<1 \Z — 

As above, we deduce by changing variables back that 

(4.8) ||5''^IUi(i/2<|f|<i) ^ ll^'^IUi- 
Also, it is trivial to see that 

(4.9) ||5'^I|l°°(1/2<|6»|<1) < £^^||t^^||L°°. 

By factoring out e in the denominator we can re- write ( [4 .71) and use Lemma 2.1 in |^ to 



deduce that: 

IT 1 — 1 f 1^1 e(n ±\jn ^ ^|| e||l/2 ii e||l/2 

~ ^ A/2<|e|<i \{zle)-ef ^ ' ^ - e "iMi/2<|e|<i)ll^? llL-(i/2<|e|<i) 

— '-^ ll"^ Wjj- 

where we have used relations ( [4 .81) and ( [4 .91 ) in the last inequality. This concludes the 

proof. □ 



Remark 4.1. It follows from ( [4.1|) that the estimate for v'^ in Theorem ^?T| is uniform in e 
and t. 

4.2. Harmonic vector fields and the cutoff function. The other estimates we require 
in order to study the asymptotic problem involve derivatives of velocity. Before we begin 
examining these estimates we must address the issue that uj'^ and m"^ are defined on an 
e-dependent domain, and if we want to explore their asymptotics, and, to use standard 
functional analysis reasoning, we must make sure they are all in the same function space. 
We will introduce a suitable ^-dependent cutoff function for a neighborhood of which, 
for reasons which will become clear in the next section, is adapted to the geometry of the 
domains. 

Let G C°°(M) a cutoff function with the properties that < < 1, is monotone 
increasing, 0(s) = 1 if s > 2 and 0(s) = is s < 1 + a, for some < a < 1. For x G IIe, set 

0^ = 0^(a;) = 0(|T^(x)|)=0(|T(x/e)|), 

and 0^(a;) = for x G fig. Clearly 0^ G C°°(]R^), vanishing on a neighborhood of ^2;.. We 
require some properties of V0'^, which we collect in the following Lemma. 

Lemma 4.1. The cutoff (p^ defined above has the following properties: 
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(a) if is the harmonic vector field for then ■ = m Ilg, 

(b) there exists a constant C > such that |V0^| <C/e, 

(c) there exists a constant C > such that the Lebesgue measure of the the support of 

is bounded by Ce^. 

In particular, || V0^||li(ir2) and || V0^||l2(]I52) is bounded as e ^ 
Proof. First we observe that 

^,=(.) = _Lv-iog|T.wi = ^-J^v-riw, 

and 

V0^ = 0'(|T^(a;)|)V|T^|(x), 

so that the first assertion follows. 
Next we compute 



|V0"(x)| < C 



^'^''^ -DT'ir) 



\THx) 



e 



by (IJ) and since DT^{x) = {l/e)DT{x/e). 

Finally, the support of is included in the set {x G | 1 + a < |T^(a;)| < 2}. The 
Lebesgue measure of this set can be estimated by: 

dx = e^dy = Ce^. 

l+a<|rE(x)l<2 >''l+a<|T(|/)|<2 

□ 

We will use the cutoff 0^ to uniformize the domains under consideration. All the functions 
and vector fields defined in are hereafter to be extended arbitrarily to the full plane (for 
instance by assigning zero value inside il^). We will only be working with these extensions 
multiplied by 0^, which makes the extension chosen irrelevant. 

We require more detailed information on the asymptotic behavior of than what was 
provided by observation (|3.7| ). This information is encoded in the following Lemma. 

Lemma 4.2. Let H = x-^/27i\x\'^ denote the basic harmonic vector field on \ {0} and 
fix R > 0. Then, 

strongly in L^{B{0; R)) as e ^ 0. 
Proof. We estimate directly: 

[ \(f)'H'-H\dx< I \(f)'H'-H\dx+ I \(j)'H'\dx 

J\x\<R J {\x\<R}n{\T<^(x)\>2} J l+a<\T^ {x)\<2 



+ / \H\dx=Ii+l2+l3- 
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We begin by noting that X3 ^ as e ^ because the function \H\ is locally integrable 
and the Lebesgue measure of the set Qe U {|T^(a;)| < 2} tends to zero, since it is contained 
in a ball centered at the origin, with vanishing radius. Next we estimate T2: 



1X0 1 < 



\H^\dx = e 



l+a<|TE(x)|<2 



\Hii{y)\dy, 



l+a<\T{y)\<2 



which clearly vanishes as e ^ 0. Above we have changed variables to y = x/e and we used 
the scaling of H"^ from ( p.6| ). 
Lastly we address Xi. We find: 



1^1 1 



{|2^l<RMn{|T(2^)|>2} 



{\y\<R/s}n{\T{y)\>2} 



\Hn{y)-H{y)\dy 



DT\y)T{y) 



2vr|r(|/)|^ 



2-K\y\ 



dy. 



Now we will use the expression for T from Lemma T{y) = f3y + h{y), with /? G M, 
P 0, h a holomorphic, bounded function whose derivative \Dh{y)\ < C/\y\'^ for some 
C > 0. We then find: 



^1 



= Ce 
< Ce 



{\y\<R/e}n{\T{y)\>2} 



iPI + Dh\yWy + h{y)y 



If^y + HyW 

iDh\y)){Py + h{y)y 



< Ce 



'{\y\<R/e}n{\T{y)\>2} 

+ Ce 
T^dy + C 

'\Tiy)\>2 \yr 

using ( pl6| ) in the second integral, 

<C{e + e\og{R/e)), 
which vanishes as e — 0, as we wished 



\Py + HyW 



{\y\<R/e}n{\T{y)\>2} 



[31 



\Py\' 

dy 

{f3y + h{y))^ 
\f3y+h{y)\^ 

\Hy)\ 



dy 



Py^ 



dy 



{|s/l<^?Mn{|r(s/)|>2} MPv + Hv) I 



dy, 



□ 



4.3. Temporal estimates. The standard way to derive temporal estimates from spa- 
tial regularity is to use the PDE directly. In our problem this is easy to do for vor- 
ticity, i.e., system and Theorem imply that, for any T > 0, urf is bounded in 



L-([0,r];W^, 



-1,1/ 



loc 



')). Indeed, 



div {{v' 



[a 



with tu^ and v"^ bounded in L°° and if^ bounded in Lj^^, see (|3.71) . However, we require 
temporal information on velocity, which is more difficult to obtain, because of the presence 
of the pressure in the velocity equation. Obtaining such estimates is precisely the subject 
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of the article [Q. The specific case of the exterior domain was not discussed in that paper, 
so we will adapt the technique used there to the present situation. 

Recall the definition of the stream function ( p.2|) and introduce the analogous family of 
stream functions given by ip^ = G'^[uj'^]. 

Proposition 4.1. For each R,T > there exists a constant C = C{R,T) > such that 



ip{x)4't{x, t)dx 



^ ^11 ,,1/2 II ||l/2 



for every in Coijle) with support contained in -8(0; R) and all < t < T . 

Proof. First we observe that there exists R> such that uj'^{-,t) has support contained in 
the ball B{0; R) for all < t < T and e > 0. To see this, first let Rq be such that B{0; Rq) 
contains the support of Uq and note that equation (E^) means that u'^ is transported 
by the velocity field u"^ = v"^ + {a — m)H^ , with uniformly bounded by a constant C 
independent of e and t and H'^ is bounded by another constant C, also independent of 
e, outside of 5(0; i?o)- Hence the support of co'^(-,t) is contained in 5(0, i?o + 2Ct), and 
taking R = Rq + 2CT will work. Additionally, uf is also compactly supported in the same 
ball. We differentiate the definition of the stream function and write: 

rt=G^[u;t]. 

This means that 

Aip^ = ul in lie, and ^/'f = on T^. 

We require information on the behavior of ipl near infinity, which can be readily obtained 
from the compactness of the support of ul and the explicit expression for G^, given by 
substituting T by T^. We obtain that, for each fixed e, 

(4.10) \il)l{x,t) - L[ujl{-,t)]{x)\ = 0(l/|x|), when oo, 

with, the functional L defined by 



for any test function C,. 

Additionally, by ( p^.8D , we have that 

(4.11) iWt I = C(l/|xp), when |x| ^ oo. 

Let (y9 be a fixed test function in Co(ne) with support contained in B{0]R). Let M 
ip{x)dx. Define 

M 

rj^G^[^] + — \og\T% 
Ztc 

Then r] satisfies the following properties: 

Arj = if in lie, = on Fg, 
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(4.12) 
and 
(4.13) 

Then we have 



M 

ri{x) = — log |r^(x)| + L\lp] + C(l/|x|), when \x\ oo, 



|Vr/| = — V(log \T'\) + 0{l/\x\^), when |x| ^ oo. 
2n 



ip{x)ipl{x, t)dx 



Ari{x)4'^{x, t)dx 



ri{x)/^ipl{x, t)dx + / i'ipt'Vr] — riVip^) ■ fids 
Jan. 



where the boundary terms include the terms at infinity. Looking at the boundary terms 
expressed in J , we see that the terms integrated on vanish, whereas from the asymptotic 
formulas ( [4.10| - [4.13| ) the terms at infinity are bounded, in such a way that we arrive at: 



(4.14) 



\J\ < CML[ijt] < Cy\\Li\L[uj 



We claim that is bounded, uniformly in < t < T and e > 0. To see that we use 

[E^) in the following way 



1 
1 

~ 2^ 



i\og\T'{y)\)u%y,t)-\/uj%y,t)dy 



\/{\og\r{y)\)-u%y,t)cu%y,t)dy 



V(log|r(l/)|) ■ (v'iy^t) + {a-m)H'iy))co%y,t)dy 



By ( p.ll|) and (|3.6| ), one way of expressing the harmonic vector field is as 



H' = —V^\og\T% 
Ztt 

so that the dangerous term in the integral above disappears, leaving: 



iH%y))^-v'{y,t)cu%y,t)dy 



< C||i/^||ii(s(o;_R))||t<;o||L°° < C, 
by Theorem [4.1| , ( |3.7] ) and ( ^.1| ). Together with ( |4.14| ), this means that 

\j\<cml^<cm];i\m];', 

since the support of ip is contained in B{0; R). We are left with estimating |X|. 
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We have: 



1X1 



ri{x)Aip^{x, t)dx 



/ Vr]{x) ■ {v^ + (a - m)H^)u^dx 



with no boundary terms because u'^ is tangent to and uo^ has compact support. Thus, 

|X| < ||Vr/|Uoo||(t;= + (a-m)i^>^|Ui(s(o;ij)) <C^||Vr/|Uoo, 

again using Theorem (|3.7|) and ( [4. 1|) . 

Finally, note that, by the definition of 77, we find that 

V^r] = K'[ip\+MH^, 

with M = J ip. Thus, we use Theorem ^TT] , with v9 in place of a;^, to conclude that: 

IIV7 II / r^ll l|l/2|| ||l/2 

II V?7||l-' < C'IIv^IIlocIIv^IIli • 
Putting together the estimates for |X| and \J'\ concludes the proof. □ 

We write the conclusion of this Proposition more explicitly in terms of temporal estimates 
for the quantities of interest in the Corollary below. Recall 0^ the cutoff function from the 
previous section and consider 0^a;^, (j)^v^ and (p^ip^ as functions (or vector fields) defined in 
the full plane. 

Corollary 4.1. Let R,T > 0. Then there exists a constant C = C{R,T) > such that 

(4.15) mV%{-,t)\\H-HBiO;R))<C, 

for all e and < t < T . 

Proof. Let C e Hl{B{0; R)) x HliB{0; R)). Then, as vf = V^^f , 



l(C,(0V),(-,t))| 



Cfvti;t) 



curl iCf Wti;t) 



curl (C)0^V^,^(-,t) + J C-VVV'[(-,t) 
< C||0^ curl a'/^U^ curl Cll^^' + C||CV0l|i/i ||CV0lli(' 
by Proposition |4.1| used first with = 0^ curl ( and then with (f = ( ■ V~^(f)'^, 

<C{\\ curl CIIloc + IICIU-), 



since || V0^||loo < C/e and || V0^||li < Ce, as proved in Lemma [4.1| , items 2 and 3, 

< CK\\h^B{0;R)), 

again by the Sobolev embedding theorem. This gives ( |4.15| ). 



□ 
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5. Passing to the limit 



In this section we will use the estimates developed in the previous section to obtain limit 
equations describing the behavior of the flow obtained in the limit e 0. 

5.1. Strong compactness in velocity. We will use a parametrized version of Tartar and 
Murat's Div-Curl Lemma to derive strong compactness in for the sequence of velocities 
(p'^v^. We will include here the precise statement of this version of the Div-Curl Lemma, 
whose proof can be found in 0. 

Lemma 5.1. Fix T > and let and {G'^{-,t)} be vector fields on M" for < 

t <T. Suppose that: 

(a) both F and ^ G weak-* m L°°([0, T]; M")) and also strongly m 
C{[Q,T]-H~l{W-W^)); 

(b) {div F^} IS precompact m C([0, T]]H~l{W)); 

(c) {curl G'^} IS precompact m C([0, T]; i/;~^(M"; A")). 
Then F' ■ G^ ^ F ■ G m V'{[Q,T] x M"). 

We will use the Div-Curl Lemma with F'^ = G^ = (p^v^. We begin by observing that 
by Theorem [4.1| we know that {(p^v^} is bounded in L°°([0,T] x M^), which is contained 
in L°^{[0,T];Ll^{M.^)). By Corollary {i(f)'v')t} is bounded in L°^{[0,T]; H'^), so that 
the sequence {cp^v^} is equicontinuous from [0, T] to H^^. Recall that Lf^^ is compactly 
embedded into H^^, so that we can use the Aubin-Lions Lemma (see [|14|) to conclude that 
{(p'^v^} is precompact in C([0, T]; H{^1). Passing to a subsequence if necessary, we conclude 
that there exists v e C([0, T]; R-l) n L~([0, T] x R^) such that 

(5.1) ^"v'^v 

strongly in C([0, T]; H^-l) and weak-* in L~([0, T];LlJ. 

Theorem 5.1. We have that cp^v'^ — > v strongly in L^q^([0,T] x R^). 

Proof. It is enough to verify the remaining hypothesis of the Div-Curl Lemma, i.e. that 
div (0^t>^) and curl {(p^v^) are precompact in C([0, T]; First note that, by Corollary 

13, both (div {(p^v^))t = div {(p'v')t 
L°°{[0,T]-H-^). Also, we have that 

div {(p'v') = ■ V(p' 



and (curl {(p^v'^))t = curl {(p'^v^)t are bounded in 



which, by Lemma ^]T] and Theorem [4.1| is bounded in L°°([0, T]; L^), and 

curl {(p^v") = ■ S/^cp" + (p^u^ 

which, by Lemma |4.1| , Theorem ^?T] and ( [4.1| ) is bounded in L°°([0, T]; L^). Since is 
compactly embedded in H^^]., once again using the Aubin-Lions Lemma, we conclude that 
both the divergence and the curl are precompact in C([0, T]; By Lemma |5]l 

infer that 



we 



10V 



in V\ 
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which in turn imphes that converges strongly to v in Lf^^, as we wished. □ 

5.2. The asymptotic vorticity equation. We begin by observing that the sequence 
{0^a;^} is bounded in L°°([0,T] x M^), so that, passing to a subsequence if necessary, we 
have that 

(5.2) (l)'uo' uj, weak * L~([0, T] x R^). 
We are already in possession of a limit velocity 

(5.3) u = V + {a — m)H, 

where v is the strong limit of (p'^v'^ from the previous subsection and H is the strong limit 
of (l)^H'^, introduced in Lemma [4.2| . The purpose of this section is to prove that uj and u 
satisfy, in an appropriate weak sense, the system: 



(5.4) 



Ut + u-Vu = 0, in X (0, oo) 

div u = and curl u = u + {a — m)6, in x [0, oo) 

\u\ — > 0, as \x\ — s> oo 

uj{x, 0) = uJo{x), in 



d2 



Above, 6 is the Dirac delta centered at the origin. 

Definition 5.1. The pair {u,uj), with u a vector field in L°°{\{]^oo)\L]^J^^)) and uj G 
L°°([0, oo); L°°(M^)) is a weak solution of system ( ^-4) if 

(a) For any test function ip G C^([0, oo) x M^) we have: 

poo p p 

(ftU! dxdt + / / V(p-uujdxdt+ / (p{x,0)uJo{x) dx = 0, 
Jo Jk2 Jk2 

(b) we have div u = and curl u = u + {a — m)6 in the sense of distributions, with 
\u\ ^0 at infinity. 

When a = m the definition above reduces to the standard definition of weak solution 
for the vorticity formulation of the incompressible 2D Euler equations. We now state and 
prove the main result in this article. 

Theorem 5.2. The pair {u,uj) given by ( |5.5| j and ( \5.3i ) is a weak solution of the system 

Proof. We will begin by verifying that {u, uj) satisfy the linear elliptic system corresponding 
to the second condition in Definition |5.1| . Recall that = f ^ + (a — m)H'^ . First observe 



that (p'^u'^ — > u strongly in Ll^^, by Lemma and Theorem |5.1| . Hence, by Lemma 
and Theorem 14. 11, 



div u = lim div {tp^u^) = \im.v^ ■ Vcff = 0, 
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where the hmits were taken in the sense of distributions. Similarly, we use that (j)^v^ — > v 
strongly in Lf^^, so that 

curl V = lim curl = \im cp'^uj'^ + ■ V"*"*^^ = co, 

in the sense of distributions. Hence, 

curl u = curl v + {a — m)curl H = u + {a — m)5. 

The velocity u satisfies the condition |m| at infinity because the convergence of 
to u is uniform outside a ball containing the origin, as can be checked directly by the 
explicit expressions for /^'^[u;^] and if^, using the uniform compact support of uj^ . 

Next, we introduce an auxiliary nonlinear functional I^. Given any test function Lp G 
C,°°([0,oo) X M2) let: 

p-oo r poo r 

Jo Jk^ Jo Jr^ 

Fix ip G C^([0, oo) X M^). The proof that {u, u) is a weak solution proceeds in two steps. 
First we will show that 

^ei^] + / ^(^^ O)co'o(a^) dx — ^ 0, 
Jk2 

when e — i> 0. The second step consists of showing that 

Ie[^]^ / / (ftuj dxdt + / / W(p ■ ULU dxdt, 
Jo Jr'^ Jo Jr'^ 

as £ — i> 0. Clearly these two steps complete the proof. 

We address the first step. As and uj^ satisfy (-Ee), it can be easily seen that 

dxdt 



JR2 



V{^{(p'f)-u'uj'dxdt- / ip{x,Q){(p'f{x)uQ{x)dx 



Thus we compute: 

isM = - 



JM2 



JR2 



e\2/ 



*^^u^)uo^ dxdt 



V9(x,O)(0'^) {x)ujq{x) dx 



-2 



JR"^ 



(fVcf)" ■ {(f)''v')u' dxdt - / (p{x,0){fy{x)uJo{x)dx, 



since ■ = 0, by Lemma 4.1. It is easily seen that 



+ 



(p{x, O)(0^) (x)u;o(x) dx 



< ho' 



V||loo||(^V0" 



0, 



as £ —> 0, by ( |4.1| ), Theorem ^]T] and Lemma [4.1| . On the other hand, since the support 
of ujq does not contain the origin, it follows that for all e sufficiently small, (0^)^a;o = uJo- 
This concludes the first step. 
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For the second step we begin by noting that the hnear term offers no difficulty. The 
nonhnear term consists of the weak-strong pair vorticity-velocity, with the vorticity 0^ct;^ 
converging in the wealc-* topology of L°°((0, oo) x M^) to uj and the localized velocity 
Vip ■ converging in the strong topology of i^H(O) ^) ^ to u. This concludes the 

second step. □ 

\i a = m the whole sequence converges to without needing to pass to a sub- 
sequence, due to the uniqueness portion of Yudovich's Theorem, see |]T^. Indeed, our 
argument shows that for any sequence cj^ ^ 0, there exists a subsequence u^- such that 
cu^'^i ^ uj and that is a weak solution of 2D Euler in the full plane with i^o as initial data. 
By Yudovich's Theorem, uj is uniquely determined, so that any accumulation point of the 
precompact sequence tu^'^ is precisely equal to u. This implies that the whole sequence uj'^'' 
converges to uj, and, as this subsequence was arbitrary, our contention follows. 

If a 7^ m. Theorem ^.2| implies existence of a weak solution for the limit equation (|5.4| ), 
which is, roughly speaking, the usual Euler equation with an embedded point vortex back- 
ground. Also, its restriction to \ {0} is the standard Euler equation (see relation ( p.8| )). 
The more physically meaningful presentation of the incompressible 2D Euler equations is 
the velocity form. Clearly, if a = m, the limit flow satisfies the usual velocity form of 
the incompressible Euler equations, but if a 7^ m it is not entirely clear what form the 
asymptotic balance of momentum equations should take. We clarify this issue in the next 
subsection. 

5.3. The asymptotic velocity equation. We have obtained a satisfactory description 
of the small obstacle limit expressed in terms of vortex dynamics. The purpose of this last 
subsection is to obtain a description of the limit behavior expressed in terms of flow velocity. 
We will avoid introducing explicitly the weak forms of the velocity equation, keeping the 
discussion less technical than before. More specifically, our purpose is to determine the 
specific form of the momentum equations in the small obstacle asymptotics, in the case 
a ^ m. The main difficulty is making sense of the term H ® H . We will present several 
equivalent forms for the limit equation (|5.4|). Recall that 



H ^ 



27r|x|2 

First we note that ( ^.4|) can be re-written as: 

{ujt + div(t>ct;) + (a — m) dw^Hoj) = 
V = K[uj] 
uj{x,0) = UJo{x) 

which is clearly well defined in the sense of distributions. This can be easily seen to be 
equivalent to ( ^.4|) since one has that div H = and curl if = 6 in the sense of D'(R^). 
Writing (p.4|) in this way will lead to a weak formulation for the velocity equation. More 
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precisely, we next show that if ( |5.5|) holds, then we have the following equation for v 

{Vt + v- Vf + (a — m) div(f ®H + H®v) — {a — m)v{{i)^5 = Vp 
divt; = 

In order to prove the equivalence of ( |5.5|) and ( |5.6| ) it is sufficient to show that 
(5.7) curl [div(i; ®H + H®v)- v{'d)^5\ = dw{H curl v) 

for all divergence free vector fields v belonging to the space VFj^f for some p > 2. Indeed, 
if ( |5.7| ) holds then we get for u = curl v 



= curl = curl ^ + f ■ Vf + (a — m) div(t' H + H v) — (a — m)t'(0)"'"5] 
= curl(f i + V ■ Vf ) + (a — m) curl [div(t' H + H v) — f (0)"'"5] 
= ujt + V ■ Vcu + (a — m) div{Huj) 



so relation ( p.5|) holds true. And vice versa, if ( ^.5|) holds then we deduce that the left 
hand side of ( p. 61 ) has zero curl so it must be a gradient. 



We now prove ( p. 71 ) under the hypothesis that divt> = and v E Wj^f, p> 2. First note 
that if curl f is well defined since curlf e Lf^^ and H G L^^^ for all q < 2. Moreover, since 
PViof C C°, f (0) is well defined, too. Next, it suffices to prove ( |5.7| ) for smooth v since we 
can pass to the limit on a sequence of smooth approximations of v that converge strongly 
in Wy^^ and C°. Now, it is trivial to che 
distributions, the following identity holds 



in Wy^^ and C°. Now, it is trivial to check that, for a 2 x 2 matrix A with coefficients 



IT A T curl Ci 
curl div A = div , ^ 
\cnn 02 

where denotes the i-th column of A. For smooth f , we now deduce that 

curl div(. 0H + H0V) = div M^""^'] + 

^ ^ \^curl(f ii2) + curl(iif2) J 

= div{H curl f + f ■ V^H + v curl H + H ■ V-^v) . 
It is a simple computation to check that 

div(f ■ V^H + H ■ V-^v) = V ■ div H + H ■ divv + curl w div if + curl H div v. 
We therefore get the following general formula 

curl div(f ® H + H v) = dw(H curl v + v curl H) + v ■ V"*" div H 

+ ii ■ V"*" div V + curl div H + curl ii div v. 
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Taking into account that div v = div H = and curl H = 5 we infer that 
curl div{v ^ H + H ^ v) = div{H curl v + v6) 

= div(iJ curl v) + div [v{0)6] 
= div{H curl v) + curl [t>(0)"'"5] . 
Relation ( ^.7|) now follows and so does the formulation (|5.6|) . 



We would now like to give a formulation in terms of u only. Simply replacing v by 
u — {a — m)H in ( p.6|) is not very enlightening. On the other hand, u is not L"^^^ since 
H ^ Lj^Qj,. Therefore, a formulation which makes use of m ® u cannot be made rigorous. 
Nevertheless, it is still desirable to obtain such a formulation in order to clarify which form 
the limit equation for the velocity takes. 

We will proceed as follows. First note that 

u 1^ u = V V + (a — m){v H + H v) + (a — niYH H. 

All these terms except H ®H are well defined. In order to give a sense to m (g) -u up to we 
will simply extend if ® -ff up to by its finite part Pf(if ® H) that we define as follows. 

Definition 5.2. Let g he a function which is homogeneous of degree —2 and of class C°° 
on M? \ {0}. The finite part of g is the following distribution on M? : 

C^iR')3^^{Ffg,^)=hm(^J g^ - m J g)- 

\x\>e e<|2:|<A 

Remark 5.1. This definition depends on the choice of the positive parameter A that we fix 
once and for all. 

We claim that if we extend u ® u to a distribution of P'(M^) by 

u ® u = V V + {a - m){v H + H v) + {a - mf Pf{H H) 
then the limit velocity u verifies the following PDE 

{Ut + div{u 0u) = — Vp + (a — m)v{0)'^5, 
divu = 
u{x, 0) = K[uJo] + (a - m)H. 

This clearly follows from ( ^.6|) if we are able to prove that dwPf{H ^ H) is a gradient. We 
now show that curl div Pf (if ® ii) = in V'(M?); this clearly implies that divPf(ii' (g) H) 
is a gradient. Let (p G C^(M^) be a test function. By definition 

(curldivPf(/f , -{Ff{H ® H) , V ® V^<^) 

= -lim( y ■ (V® V^V^) - V® V^<^(0) ■ J H(S)Hy 

\x\>e e<\x\<X 
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A simple calculation shows that / if (8) if is a matrix proportional to the identity while 

e<|a:|<A 

V (8) V~^(p{0) is a trace free matrix. Therefore 



e<\x\<X 

It is easy to check that, for x 7^ 0, 



Next, we integrate by parts and use Stokes formula to deduce that 
(curldivPf(//® if) = lim(^ y div(ii ® ii) ■ V^v? + 7 ^ {H H) ■ {x ^V^cp)^ 



|a;|>£ 



= lim(^--J y d\Y{H®H)-x^ip+^ j {H ® H) ■ (x^V^^p)^, 



\x\=e \x\=e 

where we have used that, for x ^ 0, curl div(ii H) = 0. On the right-hand side, the first 
term vanishes since div(if ® H) is proportional to x, that is orthogonal to x^. The second 
term also vanishes since 

(H(^H)-(x(^ V^<^) = (H-x)(H- V^93) = -^T^(x^ ■x)(H- V^^) = 0. 

27r\xf 

This concludes the proof of the relation curldivPf(if (g) if) = 0. Finally, we observe in 
passing that 

div Pf (ii (g) ii) = ^ V (Pf ^ + nS). 



\x\ 



6. Final remarks and conclusion 

The results we have proved here are not very surprising after the nature of the harmonic 
part of the flow has been clarifled. In fact, a good illustration of the results obtained can 
be explicitly computed if one considers the limit flow associated to the motion of a single 
point vortex in the exterior of a disk of vanishing radius. To be precise, consider the motion 
of a single point vortex of strength m > in the exterior of the disk B{0; r), r < 1, initially 
located at (1,0). Let Pr = Prif) denote the trajectory of this point vortex, which will 
remain on the circle of radius 1 around the origin, moving with constant angular velocity. 
At each flxed time the velocity fleld can be computed using the method of images as: 

m {x — Pr(t))^ m {x — r'^Prit))-^ a x^ 



Ur(x,t) = — — — - TT-T— + 



2n \x-Pr{t)\^ 2tt \x - r^Pr{t)\^ 27i\x\'^' 

This flow is equivalent to flow in the full plane associated with three point vortices: the 
original one at Pr{t), the image point vortex at r^Pr{t), with strength — m and a point 
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vortex of arbitrary strength a at the center of disk. The location of the image vortex is 
the inversion, with respect to the disk of radius r, of the location of the original vortex. 
The point vortex of arbitrary strength at the center is associated with the choice of the 
harmonic part of the flow. Curiously, when the method of images for flow in the exterior 
of a disk was discussed in |T^, Saffman simply assumed that a = m was the correct choice 
of harmonic part. In the case a = m, we have that 

so it is easy to verify that the angular velocity of Pr{t) vanishes as the radius of the disk 
vanishes, so that Pr{t) converges pointwise in time to Pr{0) = (1, 0). Furthermore it can be 
readily checked that the velocity fields Ur converge to the velocity field due to a single point 
vortex in full plane flow as r — if and only if a = m. In short, what we accomplished 
in this paper is to verify that the evanescent obstacle exterior flow, for a general domain 
geometry and general vorticity, behaves exactly like the corresponding limit of point vortex 
flow in the exterior of a vanishing disk. 



Fig. 1 Curves Xr{T., T) for T ^ 2 and r ^ 0; 0.025; 0.05; 0.1. 



One interesting feature of the limit process in the case a = mis that, although both the 
approximating and the limit flows are smooth, the particle trajectories do not converge 
uniformly. To see that, we use the illustration described above. Set m = 1 = a and let 
Xr = Xr{Q,t) denote the particle trajectory under the flow Ur, starting at the Lagrangian 
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marker Q, and let Xq = Xo(Q,t) denote the particle trajectory under the limit flow. Let 
S = {|x| = 1/4} and assume r < 1/4. Clearly, for any t > 0, and Xo(S,t) are 

smooth Jordan curves in the plane. It can be checked that there exists T > such that 
the origin lies outside Xo(S,T). Fix such a time T > 0. On one hand one expects to 
have Xr(E,T) converging to Xo(E,T) as r — »• 0. On the other hand, the origin lies in 
the interior of Xr(S,T) for all r > 0, as the flow has to remain in the exterior of 5(0; r), 
see Figure 1. This means that the convergence as r ^ of the maps Q ^ Xr{Q,T) to 
the limit map Xo{Q,T) cannot be uniform, even in the compact parts of their common 
domain. 

Below we add two additional remarks concerning the results we have obtained. 

(a) In the case 7 = a — m 7^ 0, one may ask why fix the circulation 7 around the small 
obstacle, since this implies very large, maybe unphysical velocities at the boundary 
of the obstacle, and not consider it some appropriate function of e, perhaps vanishing 
when £ — >■ 0. In fact, 7 plays the role of a passive parameter in the argument we 
have presented, so there would be no change in the argument if we consider 7 to be a 
function of e. The limit flow would depend on 7(5) only through its limit when 5 — > 
in precisely the same manner as presented. 

(b) For the sake of simplicity, we have presented our argument for smooth, compactly 
supported initial vorticities, but the argument can be easily performed for compactly 
supported vorticities in L^, p > 2. The argument does not work for p < 2 because 
there are serious difficulties in making sense of the term Huj in that case, and the 
value f (0), which appears in the limit velocity ( p.8| ), also becomes ambiguous. The 
case p < 2 is thus an interesting open problem. 

Let us point out some of the natural questions raised by the research presented here. 
First, an analogous question can be asked with regard to the 2D incompressible Navier- 
Stokes equations, and this is the subject of current investigation by the authors. For 
incompressible 3D Euler, the problem initially becomes proving that a smooth solution 
exists for a time that is independent of the size of the domain, something we did not 
investigate, but that appears to be difficult. Of course, the most interesting situation is 
the same limit for 3D Navier-Stokes, but it makes sense to work out the 2D equations first. 
Once the viscous problem has been understood, one may ask about the interaction of the 
small obstacle and small viscosity limits, which in two dimensions is a simplification of 
the classical open problem of convergence of Navier-Stokes solutions to Euler solutions in 



the presence of boundaries, see |]15| for an account of the state-of-the-art of this problem. 
Another natural question is to understand the same limit with more than one obstacle, 
with perhaps only some of the obstacles vanishing. Finally, it would be interesting to 
obtain a description of the leading order correction associated to the evanescent obstacle, 
in the case a = m, with respect to the unperturbed full plane problem. 
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